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CHAPTER I 


NTRODUCTION AND SUMMARY 


1. Introduction 

The steady state spatial distribution of neutron fluxes in matter is a 
problem of prime importance to the designer of chain reacting piles, and 
has been the subject of very much theoretical and experimental work. The 
case in which the neutron energy is large compared to that of the nuclei 
constituting the mediun (the "slowing down" case) is of first importance 
and is rather extensively investigated. The alternate case in which the 
neutrons are in energy equilibrium with the medium is of lesser importance 
and comparatively little experimental work has been done, although the 
fundamental theory is well in hand. The cross-sections for neutron scat- 
tering and absorption appropriate to the first case are the "ordinary" 
energy-dependent cross-sections (or simple functions of these) which have 
been so extensively investigated by various neutron spectrometry groups. 
The cross~sections appropriate to the second case are more vaguely defined. 
Perhaps as adequate a definition of these parameters as any is the "ad hoc" 
one - those numbers which, when inserted in the mathematical theory of 
particle transport yield proper solutions for the neutron flux. It can 
be demonstrated that these numbers are not in general the average values 
of the energy dependent cross sections weighted by the appropriate energy 
distribution function for the ambient temperature. 

The "ad hoc” definition suggests a means of experimentally measuring 
the numbers. One requires a medium - of simple geometry so that the trans- 


port problem is soluble - in which a flux of thermal equilibrium neutrons 








is maintained and in which the flux can be measured at various points with 
neutron detectors. Since two parameters are to be determined, two independent 
solutions for the flux must be available - but this is easily provided for 

by considering (1) a region "remote" from discontinuities in which the 
solution is asymptotic and (2) a region adjacent to a discontinuity - say 

a "black" boundary. One measures the flux at a sufficient number of points 

in each region to define the spatial functions. The transport problem is 
solved many times for various values of the parameters, and if fortune is 

good one pair of values will be found to yield solutions which fit the ex- 
perimental functions in both regions. 

To particularize = the chosen medium is normal water, the geometry 
one-dimensional consisting of a tank many mean free paths in its dimensions, 
traversed by a wniform thermal neutron current arriving at one face and 
absorbed by a layer of cadmium at the opposite face. The neutron flux is 
measured at points along a line normal to the cadmium facing and extending 
from the facing well into region (1). 

The one-dimensional transport problem is integrable in closed form for 
region (1). This solution is matched to the experimental data, yielding a 
characteristic equation in say the total cross-section and the scattering 
cross-section - thus the scattering section may be eliminated. The trans- 
port problem is then integrated numerically in region (2) for various 
values of total cross section and that value selected which yields best fit 
with experimental data. CAtereanmeely, this straight-forward procedure 
could not be carried out. It was found that numerical integration produced 
a function which fit the experimental data astonishingly well, but that the 


values of the parameter which yielded this good fit were variable with 








respect to the integrating procedure. Further, it was found that each 
improvement in the integrating plan led monotonously to a reduction in 
the total cross section, indicating that the process of improving the 
operation was not near convergence. The rapid increase in computing labor 
made it necessary to abandon this process before a firm solution was ob- 
tained. 

The Wiener-Hopf method, treated by Frankel and Goldberg (3), offers 
a general solution to the one-dimensional transport problem in a single 
semi-infinite mediun. However, it is available only as an exceedingly 
difficult double numerical quadrature which could not be handled for lack 
of computing facility. A solution is given for the flux at the black 
boundary only, requiring only one numerical quadrature. By equating this 
solution in terms of the paraneters to the experimental value of the 
"contact" flux, unique values are obtained for the parameters. These 
values extrapolate the trend of the numerical solutions, and are in that 
sense to be regarded as supported by the numerical method. 

The convergence on the Wiener-Hopf solution may be illustrated by 


the successive values for the macroscopic “total” cross section. They 


ares Numerical - 2.9 om7!; 209 a and 2.93 om); Wiener-Hopf - 2.10 em7! 








2. Summary 


Findings may be tabulated as follows: 














Parameter ~- Water at 24° C. Values 
Diffusion Length 2.67 + .02 cm. 
Transport Mean Free Path 0.48 + .01 cm. 
Absorption Mean Path (Equivalent) 44.2 1.3 cm. 


Microscopic Transport Cross Section 62.4 + 1.5 barns/mol. 


Microscopic Absorption Cross 


Section (Equivalent) 0.68 + 0.01 barns/mol. 


Ratio - Boundary Flux + Asymptotic 


(Exponential) Flux 0.186 
Ratio - Boundary Flux + Asymptotic 

(Hyperbolic) Flux 0.87 
Extrapolation Distance - Xo 0.32 cm. 
Ratio - Extrapolation Distance/Total 

Mean Path 0.67 


| Comparison of the cross-section with the energy-dependent values for 
protons in water by Bacher(1) and Melkonian (7) indicate an "average" 
neutron energy for transport of perhaps 0.05 electron volts in contrast 
with the Maxwell-Boltzmann mode of 0.025 e.v. This is not surprising since 
it is expected that transport is largely effected by the high energy neutrons 
which have longer free paths. In addition it must be remembered that the 
neutron flux is weighted by a velocity term so that the average flux corres- 
ponds to a higher energy than does the average neutron, and transport weights 
the cross-section by the flux rather than by the neutron density. The 1/v 
absorption in hydrogen provides a preferential sink for low energy neutrons, 
thus the thermal distribution is displaced from the Maxwell-Boltzmann toward 


higher energy. 


CHAPTER II 


TRANSPORT THEORY AND WETHODS OF ANALYSIS 


To indicate the assumptions involved and to portray the physical model 
which was analyzed, it may be well to review the fundamental theory of 
neutron transport. It must be emphasized that the model of neutron trans- 
port for neutrons in energy equilibrium with the medium, at energies well 
below the atomic binding, differs radically from that for energetic neutrons 
(E, > 1 ev) in media at normal temperature. 

In the latter case the energy transfer in neutron-nucleus interactions 
is uniformly from neutron to nucleus, the nuclei interact as free particles, 
and since the atomic vibrational energy is comparatively small, the nuclei 
may be regarded as initially at rest. None of these statements hold for 
the former case. 

In the present problem, by definition, the net average energy transfer 
between neutrons and medium is very small and due solely to the preferential 
absorption of the slower neutrons. The protons are bound in the molecule 
with energy larger than the neutron energy and the nuclei cannot interact 
as free particles. The thermal motion of the atoms and of the systems in 
which they are bound is not negligible, and so affects not only the mechanics 
of the collisions but also the frequency of occurrence since the collision 
probability is a function of the flux of nuclei as well as of neutrons. 

Inelastic interastion of neutrons with rotational and vibrational 
levels of the water molecule are probably present. The rotational levels 
are quite numerous and begin at wave number 23.8 em7!, (about .003 electron 
volts) according to Wu (9). 


Quantitive discussion of the scattering mechanism in a system as 








complex as this is not in the scope of this paper. However to avoid mis- 
interpretation of the results, it is necessary. to consider at least the 
question of whether the scattering is symmetric, that is whether all values 
of the angle of deflection of the neutron in the Laboratory System are 
equally probable, when the neutrons are in energy equilibrium with the 
mediume 

Referring to Fig. (1), 4% and ¥ are the velocities of neutron and 
molecule respectively before collision, in the laboratory system, and Ww 
the angle between them. To transfer to the center of mass system a vector 
velocity V is added of such magnitude and direction as to reduce the total 


momentum to zero. Using unity as neutron mass and M as molecule mass we 








have: 

Uae SAE as M(v cos w + Vy) = 0 
aad Vy + M(v sinw + W) = 0 
saad Vy = oe vycosw — = va 
and <a M “ 

Vy ase oe 


where v and Vy are the scalar components of VY. 


The transformation vector is thus seen to resolve into two vectors, one of 


M 
M+] 
nitude __| Ve and opposite to ¥% in direction. The vectors Vv,’ 


M+! 


and V’ are the velocities before collision in \.:. syste: and W%” and v” 





magnitude v and opposite to VW in direction, and one of rag- 





are the C.M. velocities after collision. If the scatterin 5 spherically 
symmetric in C. ie system, all values of cos YW, the deflection angle, are 
‘shahle, ico obtair tne neutron velocity YU, witer collision, in 


A a7 
CQtas AY 


6 





Laboratory system, we add (-V) vectorially to %“. The probability of 


w is sinw, and is therefore symmetric, so on the average the contribu- 





tion of ve V involves no preferred direction. However, the contribu- 
tion of oe : ¥, inserts the direction of the incoming neutron as a 
+ 


component of u . Thus the scattered neutrons will emerge from collisions 
preferentially in the direction of their pre-collision motion with a weight 
increasing with neutron energy. 

The preferentially forward scattering, as indicated, correlates rather 
poorly the direction of each path with the preceding one, less well with 
the next but one preceding, and after a large number of intervening scat- 
terings the correlation vanishes. One may attempt to obtain the number of 
neutrons passing through a small volume increment (dV) by summing up those 
neutrons scattered in all other volume increments (ae in the space which 
reach the first increment in free flight. However, since the number of 
scattered neutrons leaving dV ‘in the direction of dV depends on the aver- 
age direction of the neutrons entering dv’, this average direction must be 
found by another summing up of previously scattered neutrons and so on. 

The mathematical analysis of such a situation yields infinite order multiple 
integrals, so simplification is necessary. The difficulty does not occur if 
the scattering is isotropic, so an obvious suggestion might be to approxi- 
mate the system with a model in which the scattering is isotropic. 

If one considers a neutron scattered isotropically at a point in a non- 
absorbing medium and asks for the probability of finding it at any point in 
the space after the lapse of a long time, it is clear that the resulting 
probability field is symmetric about the initial point. Analogously, one 
considers a neutron with initial motion ¥% scattered at a point anisotrop- 


ically in the Lab system but isotropically in C.M. system. 


The probability of finding the once scattered neutron at position Ff 


is then symmetric about the moving center of mass whose velocity is 


| 
M+| 


neutron at fF is only azimuthally symretric since each spherical shell of 


Vo - However, the probability of finding the second scattered 


second scattered neutron probability is immediately again retarded by the 


factor 4 » corresponding to the motion of the center of mass of the 





second scattering system. Thus the second scattered neutrons which have 
made long flight after first scattering will be found on the average farther 
along in the direction of Vp than those which make short flights. In 
spite of the forbidding complexity of the picture some information can be 
obtained about the probability distribution after a large number of collisions. 
The path lengths z » the azimuthal scattering angles al , and the cosines 
of the polar scattering angles p » are every one statistically independent, 
and may be averaged independently. The centroid of the distribution is 
obviously found in the direction of w% and may be located by first averaging 
over the azimuthal angles starting with the last flight. One obtains 

hn 
| ree oP Ay Pap +... Xf, I], pote) . Averaging over Ae and 


Pn one has 


f = ff pP prt Pop.) Apa, ad, ff fu Pou) Mz Ph.) ue Pu.) du,du, di, +... 


The required first moment is then ), (zx + ja’ ae A oe pee . ) 


which is the expansion of » One concludes then that the average 


| maw = 
position, after very many scatterings, of a neutron projected from a point 
I= 


initial point. This parameter is called the transport mean free path. The 


from the 





in a medium is at a distance A,(I+ ie? or 


higher moments of the probability function are not available in general. 








F,. Joliot (5) has computed the probability function itself for the nth 
isotropic scattering in a quite complicated series of Hermite polynomials. 
One may now reasonably postulate that the scattering mechanism in 
an anistropically scattering medium may be rather well reproduced by an 
isotropic model, provided that the scattering length be weighted by the 
factor 1 . This ignores the defect that the higher moments are not 

ee sc ueae 

It must be emphasized that capture has not been considered and that 
the model is not correct for large capture. It then remains to be demon- 
strated that the model is satisfactory if capture is small but not negligible. 

The simplified transport model now may be more clearly defined. It 
consists of a three dimensional array of fixed lattice points at which the 
neutrons scatter isotropically, or are absorbed. 

The scattering length in the model is the transport free path in the 
real medium and is in general not known. The absorption free path in the 
model is not in general the same as the absorption length in the medium but 
is probably nearly the same if absorption is small. A more useful pair 
of parameters are the inverses of the free paths - the macroscopic cross 
sections - defined as the probability per centimeter of path of the defined 


event. Io prove the relation: 


~ap= [ 2 dz Where I is the beam intensity 
d= 2 e722 
Let i = | ; ae = oe 


The mean value of 2 is: [s eth df 
A= 


ch aly 
joer df 6 





A macroscopic absorption cross-sectior is similarly defined. 


The following notation is to be used: 


n - Neutron density, om75 


Mean macroscopic transport cross 


c+ 
1 


p 


Mean total cross section, om71 


- Microscopic total cross section, 


- Distance, cm. 


wm RG zmMyMnms 


- Distance in mean paths 


Position vector 


ss! 
1 


- Mean macroscopic absorption cross section cm 


- Neutron scalar flux density, called simply "flux," cm7* sec 


section, oma) 


-1 


~ Lattice point density of the medium, om7% 


one 


One asks for the number of neutrons scattered in incremental volume 


at position F and crossing unit normal area at 
scattered at F is YrFyL,dv - all cross 
of fF in the homogeneous medium = and of these 
about solid angle w. Since the area subtended 
Cr’ = r)* dw » 8 division by this factor is 


beam in any dw is attenuated by absorption and 


r’ per second. The total 
sections being independent 
eal ee leave fr indw 


4n 


at * 


¢ 


by dw at fF is 
indicated, and since the 


re-scattering as 


-/ = as, =_— 
ltl | ome [7 - | multiplication by this last factor 


9 


is required. The required number is then Wr) Ze ¢ ss ae = 7 | dv 


4n(¥’-7) 


The neutron flux at fF’ is then obtained by summing up the neutrons 


scattered at all fF and crossing unit area at r’ 


Sele 7 
fs Dir) € 
PF) zt Wf (F' =P 


. Thus 


dv 


-1 


In the protlem at hand, the medium is considered to be a semi-infinite 
slab - infinite in horizontal extent and downward, its upper face bounded 
by a plane from which no neutrons return, and having a finite neutron current 
passing upward through it. The coordinate x is measured from the "black" 
plane and normal to it. In this system the flux is a function of 2% alone. 
One requires the number of neutrons scattered in an incremental slab 
dz at x and passing in free flight through unit area about a point at x! 


Referring to Fig. 2, this quantity called Kiix'— xd is found as follows: 


e x) 2 
K dx a 





(2x rdrdx) 





a St (x’— on = 


rdr = pdp 


i il 
3 -ip 
K = Pots > dE) 
z |x'- | 


ll 





The integral is termed Ee, (Z|x'=xI) » the exponential integral, 
and is identical with the function —E; (x) as tabulated by Jahnke & Ehmde 
and others. 

The flux at x' is obtained by summing the contributions of all increments 
dxat x as follows: 

Pe) = = if Pre) Er (ep—x a 
Space 
This linear integral equation is perhaps more conveniently handled in mean 


free path dimensions as follows: 


ite ees 
Ye’) = Af Yee) Es Geran 1 


For the configuration space described, the limits of integration are zero 
and plus infinity. However as -~ = —~ 15 the contribution to the 
integral between negative infinity and zero becomes negligibly small, and 
the limits may be taken, with vanishing error, over all values of Z2 . 
In this asymptotic form the equation has a simple hyperbolic solution. 

The usual method of treatment, to indicate the form of the asymptotic 
solution (as in Frankel & Nelson (4)), involves transforming the origin to 
write: 

7 oo 
Ge) = e-/ Vere Fed? 
—00 


and expand Y (2'+2z) in Taylor's series about Zz’: 


; = 2 ms 7) 2 “ 
Pe’) = a 3, (121) [ Ore) tet D (2) + olde 


=e 


co 
Now = m is the mth moment of the kernel 
m, = Z E, (21) dz 
= 00 


We 





about its origin, and this vanishes by symmetry for m odd. One has then 
= ot ‘ mM “ 
Ye) = 2 ipl P (z’) yo ve" 


neglecting terms above the third order. Rearranging, one has the one- 
dimensional differential equation of elementary diffusion theory - the 


Helmholtz equations: 


" o ny = 
e+ (a mE)Pr° 


Ze 





for which the solutions are hyperbolic if TN, < » Circulars 
‘s 
mM, > 2 » and linmeer ait mM, = 22 » Now 


Zt 2 
60 
mM, = i Evin 92 = 7 , and if absorption is present, then LE >Z4 , 
so a < az - Thus we expect hyperbolic solutions, so a 
Mo =2 
trial function Pe) = a is substituted in the integral equation, and 


4 
the integration performed at the point # = 9 


ee fue 
| Diggs * Evqzeanece 


= { f rRte da + fick E,(e) de} 


Zt {fe eae + [oh Eve | 


eo 
= 2- [cosh (RZ) E,(z) dz which is integrated 
° 


by parts after putting Sy (2) in integral for. 


! = [cosh (atz)ae [—< ay 
fo] Zz 5 


13 











si Zt a =z 
Ev) | ot BE [ sinh tay dz 


The first integral vanishes for R¢€} since at upper limit E (2) is of 
-Z . 
order & and at lower limit Sint KE goes as # and the singularity of 


Ey (z) is logarithmic. Writing the second integral in exponential form 


te Ee fer BURN g-2UEX) yy 
ZR . z 


and separating the range of integration into intervals: 


Kr Mer oo © 
Pike z(I-2) eB ita) -Z(I-#) -2(I+2) 
i= te || “2 * —[ ites ae ak eo 
(< past) ! 
i-K 


Tt 


the third and fourth integrals are seen to be each E,(1) and cancel, thus 


|= pete 9. +4 





ft (-a# 
oe a = Wes 
2 7 An ae tanh"! # 


expanding in power series 








é 
oan 3 5 


2 
Inverting the series in R” we obtain 


, ep 
Oe 2 (2) Sees () ao 


Note that 2 is an eigen function in rean free path measure and the 





~KX 
function in cm. measure is © where K= R x : 


14 





Given an experimertal value of K (which is simply the logarithmic 
slope of the flux distribution plotted against X ) one may choose a value 
k 
of 2 » compute # by k= 4%, and using the series in R? (or any conveniently 


obtained solution) compute ave » Returning to the integral equation 


ee ee oe ; 
Ye) = oe f Pee’) Evga ay 42 


one may insert the computed value of 246 - It is now required to com- 
pute the eigenfunction Ye) corresponding to the selected value of 2 . 
A choice of method is available. 

LeCaine (6) has treated the problem by a variational method, presenting 
tabulated values for the flux normalized to unit flux at the black plane 
for various values of = - A rather rough solution for the parameters 
is available by using sacs eee and indicates that 2 lies between 
2.0 om~) and 2.4 cm), 

It was then decided to attempt solution by a direct numerical method, 
which was carried out to completion. In theory the method goes as follows: 


Given the integral equation as before 


Ya) a Fe f Piz’) Es igi-zi) 9@ 


The integral is replaced by a Riemann sum equation 


V3.) = Braz) Fen Es Oz > 2,)) 


which may be put in matrix notation 


15 





: z+ 
Wl 2e0 > |= 2 ee 


where ews is the absolute value of the exponential integral of argument 
(i-n) Zoe « At the argument i =n the exponential integral is singular, 


so a mean value was computed for this range. One requires: 


o 2 SF 7 _y? 
2 2 ee ive t 
az J, Evzy9@ = 2 [be 3 ay 

oo = eo - vA? 
. FEN 7 7. ea i ae 
aff: dz= ef ye Jay 


(1 ~ A) + E, (42) 


The integrating scheme is the equivalent of the trapezoidal rule but 
employing the function of mean argument instead of the chordal mean as 
the summation element. 

One proceeds by inserting a tabular trial function, perforring the 
matrix operation as indicated and emerging with a new and bettered function 


which is iterated until the variation produced by iteration is insignificant. 


16 





In practice minor departures from the straight forward process were intro- 
duced to reduce the labor. The principal one of these consisted in iterating 
for only about 10 terms of the flux vector and writing the complete vector 
by smoothing the first differences. 

A variant of this scheme was considered, which involved writing the 


equation in a perturbation forms: 


x ao 
Pata) 52) = 3 lEs(az1)~ Eua-apll Pa(z) ~ 5 (@)| 2 


62 


where Da, is the asymptotic solution, § is the boundary depression, and € 
is the perturbation of the operator E, required to extend the integration 


over all Z . 
oe 


[Ever 7 Ee’-a))| Siz) TE ¢ dz 





Va eee 
2) “Oy 

—00 
or in Riemann sum form 


= 2+ : 
Ly Wa ae) EJ iaj-2, 5 * ECan) P 


The equivalent matrices, all infinite order are: 





SBSSiis 
is) 
SBSiis 
cs) 
i 
mm 
AN 
Ad 
a 
m™m 
2 
mm 
ws 
» 
I 


zhi. §= €= - Ena ae Ey irs 
Q; 0, 5, =- Bey a 0 0 0 
Pa b, i: Ee oe. 0 0 0 

2 / = ~~ 09:0 20 
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050.0 OS OOM Os 
0, 0) 020 0 0 Omer 
O 0 00. 0) 20 ssCraG 
0 0 0 0 C50 S0mG 
E- iJ = | 00 0 Ew Ey Es my = 
0 0 0 Ea ay en re 
0 0 0 bf: 58.) Eee 
000----- 
0 0f0°— =a 


The first method has the disadvantage that the term products converge 
rather slowly for increasing Z so that each term of the iterated vector 
represents very many term products and the computation is quite laborious. 
The second method avoids this trouble since § vanishes strongly with in- 
creasing 2 , but § varies so rapidly near (2=0) that trapezoidal inte- 
gration is valid only for a very small interval Aw . On balance, method 
one appeared most feasible and was chosen. A value of 2.9 noe was obtained 
for 2 with a very fine fit of the solution to the experimental data. 

To investigate the validity of the integration, the process was then 
repeated with an improved kernel. In this case the matrix elements were 
obtained by computing mean values for the exponential integral with a unit 
weight function for each Az,. This process amounts to correcting for 
variation of the kernel over the interval but ignoring the variation of the 
function. A value of L = 2.5 om! resulted. 

A further improvement in the kernel was then tried. One assumes that 


the function ray be linearized with small error over the interval A@ . 


18 





The matrix equation is then writtens 


z ae n 
Pe.) ae fl ve AE Ege ae ey 


where Z, is the lower boundary of the nth interval, AZ is the width of 

the interval, A n is the difference [ = ] and 
; P ( rei) Piz.) . 

e is the difference [z -Z,| - The equation may be separated into 


two terms as follows: 


5, Az 
Qe.) = $e) [Pen | Ei (lz,-zn +61) °° 


Az 


An 
vara | Evia -ari et | 


The first integral (letting y = [Z; -Z,+€] = |K,+€]) 


Kn+t 
= -K —_ —Knes — 
tn oft Evy Ay = 6 n & ine EK.) + Ke Ei (K,,,,) 


The second integral, dropping absolute values: 


Kn 
-[ (y- Kn)Ei(y) dy 
K 


n 


E [etek i C(t Ka) FRE tenag Me Ee) |eeme 


19 constitutes the terms of an operator which operates on the function Yn 


and I, an operator which operates on the relative increment of the function, 


me 


19 


The process yielded again an excellent fit to thedata but for a 
value L = 2.3 cm™!. Since the computing had become impractically tedious, 
further work on this line was dropped. 


Frankel and Goldberg (3) list a solution for the boundary flux as 
: $e - |) 
Yo) = 2 \Epi- x7] 


where the asymptotic flux goes as the hyperbolic sine vanishing at an 


"extrapolated" boundary thus: 


Ya(z) = sinh (2 +Z,) 


The extrapolation distance is given by: 





J 
] a w/) 
2.0 flex tor" (gpa ge a 
“2 ~[ \—x*s* tanh’ 5s LAE, 


The characteristic equation of the system relates KR to An » so the 


boundary flux may be expressed as a function of & alone: 


eae ee zt me 
2 Yoo) 4 L(1-&) | 





=—(|-X’) = 

















] 
zy +2. 
2 Rr x4 ( 
- BS +e tee) = —— 
ses ° ) I+ 2% 
We. oka 2" | ee 
3 ie (Smee SS Ta Hyanl I+ 2D) 
2 
2. x" Nee fee = 5% 
A Se a ag ee 1 +2, 





The procedure is as follows: 


l. 


Ze 


The extrapolation distance is computed by Simpson's Rule for 
a few values of & in the neighborhood of 0.15 = 0.20 as in- 
dicated by the numerical solution, and plotted against Xx. 
Selecting a value of & and the corresponding Zo , the 
experimental boundary flux is normalized to go as 

sinh & (z + Z.) at large argument. This new normali- 
zation is obtained by halving that flux which was normalized 
to go asymptotically as Paik and multiplying by e*%o : 
The result is Ye) . 

Insert the Yo) in the series equation and solve by trial 
for mae ° 


Using the computed values of XR » return to stage (2) and 


iterate the process. Convergence is rapid. 


The described process yielded XZ = 2.10 2 0.05 om7l. 
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CHAPTER III 


EXPER IVENTAL PROCEDURE 


The medium in which the experiment was conducted was the approxi- 
mately 6' x 6' x 6' water tank which is mounted on the thermal column of 
the Oak Ridge Graphite Reactor. A current of thermal neutrons passes 
upward through the tank after filtering through the graphite column below. 
The water content is ion-exchange demineralized, and its nuclear properties 
are substantially those of distilled water. 

A disk of aluminum, faced with .020 inch cadmium, and approximately 
10" in diameter, carried on aluminum legs about 10" long, was lowered with 
a rod to a reproducible central location on the tank bottom. A foil holder 
of aluminum and lucite, fastened to the cadmium face was provided with 
adjustable mountings for positioning foils along the central normal to the 
plane. A foil was also cemented to the center of the cadmium plane for 
each exposure. The foil mountings were of lucite of minimum section cone 
sistent with stiffness, and a film of paraffin was employed as an adhesive 
to mount the foils, thus introducing no non=hydrogenous material in the 
immediate vicinity of the foils. The equipment is shown as Fig. (3). 

Indium foils were used, weighing 25 milli-grams nominal, cut from 
50 mg/om? stock. These were cleaned of the adhesive wax after each use, 
weighed and allowed to decay for at least 48 hours before re-use. 

Foil counting was done with an end mica-window Geiger-Mueller counter, 
connected through a plate-loaded pre-amplifier stage (of gain — 12) toa 
scale of 64. Counter and pre-amplifier were enclosed in a 2" lead shield, 
which reduced background to a negligible fraction of the foil activity. 


Sufficient counts were obtained in each counting run to reduce statistical 
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fluctuation to a negligible level. A number of determinations of counter 
resolving time were made by following the decay of 54 min indium 116 over 
about 5 hours. An average value of 240 psec. was obtained by solving 


the following: 


A, = Az 





| —-N,? }-N,7 
where A = extrapolated activity 
N = mean counting rate, sect 


7 = resolving time, sec. 

A correction factor L/l-N> was plotted against counting rate and 
applied to all counting runs to correct for coincidence losses. 

For exposure, the foil holders (after mounting the foils) were adjusted 
approximately to position with respect to the cadmium plane using a common 
scale. A vernier height gauge was then used to determine the exact foil 
position (estimated probable error * .05 mn.). The plane and foil holder 
assembly was then lowered quickly to the tank bottom and exposed for times 
ranging from 5 to 15 minutes depending on foil position to obtain reasonable 
counting rates. No attempt was made to measure either exposure or foil 
"cooling" times accurately since one foil always and two foils usually were 
exposed at fixed positions to monitor the activity of the variable foil. 

The positions at x = 0 and x = 5.18 em. are the monitor positions and data 
at these points represents an average of about 15 measurements. In addition 
to the exposures made with axial foils, two exposures were made with foils 
set in a plane parallel to the cadmium to demonstrate that the geometry was 
actually one dimensional. A plot of the diametral variation of flux appears 


as Fig. (6). The maximum deviation from the axial flux amounts to about 


235 





5% reduction at 8 cm. radius. This is considered satisfactory. To test 
the flux for fast components, a cadmium difference on Indium was taken. 
The ratio of activity in Cd. to activity in the "shadow" of the cadmium 
envelope (a most conservative condition) exceeded 5,000 - therefore foil 
activation due to the 1.44 volt resonance may be safely neglected. 

To analyze the data by the numerical method, it was first plotted 
against the coordinate x in cms. on semi-log paper and the fair curve 
drawn by eye. A clearly exponential branch was obtained beginning at 


x~4.8 cm. The decay constant k of the exponential was solved by 


Ye) 
y Sn (p (b) 
a —"b 


and fixed at k = 0.375. The scale of the values of (D(x) was then adjusted 


oe A trial value of 2 was then selected, 


so that (D (asymptotic) = 
and the data plotted, faired and tabulated as a function of # in mean path 
measure. The eigenvalue = was computed, the Es (mA!) tabulated 
on a strip for Az e2, and numerical iteration of the equation carried 
out until it was obvious that the solution near the plane at # = 0 would 
not fit the data. The trial was then terminated and another trial made for 
a different value of 2 . 

As the solutions were improved it became clear that a correction was 
necessary for the remainder contributed by the terms beyond the last tabular 
value of the BE,. The correction was inserted by computing the remainder 
with respect to a proper exponential function and dividing by the value of 
this function at the first untabulated term of the FE}. This number is 
then used to terminate the FE, tabulation. Since the last tabulated value 


of the E; was for mAZ = 7.4 one obtains: 
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rad 
a7.68 
é ee 


7S 


which may be integrated by writing the E, as an integral, then inverting 
the order of integration to yield 

seiet (pure ee 

z © {-€, (75) * RS Peat 
which must be computed individually for each value of BAe 

Iteration was carried out by computing only every fifth term of the 
trial function, then interpolating the new function by making the first 
differences regular. The final integration was done in detail in the 
region close to the black boundary. 

Plots of the flux as a function of distance from the black boundary 
appear as Fig. ‘4) and (5). The curves were obtained by using the Le Caine 
results to interpolate the data near the boundary, and the numerically 
iterated function in the region more than 0.5 cm. from the boundary. The 
latter process is considered reasonable since the solutions in this region 
were invariant to changes in the integrating scheme, they are constrained 
by the approach to the asymptotic function, and they fit the data quite well. 

Desirable extensions to the analysis might include: 

(1) Computation of the Wiener-Hopf solution for perhaps 15 or 
20 points to firm up confidence in the fit. 

(2) Extensive analysis of the numerical iteration method to estab- 
lish the level of computing detail necessary to yield accurate 
solutions using the true values of the parameters. 

Both propositions require computing machine service to reduce program time 


to a practical value. 
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APPENDIX II 


The flux is tabulated as a function of Z for two values of the 
total cross section. These are solutions corresponding (for -Es 2.3) 
to the most elaborate numerical treatment, and (for 2 = 2.1) to the best 
reconciliation of Wiener-Hopf, variational, and numerical and experi- 
mental solutions. 


Zo E221) Fe 2.5) ZW Se 2.1) d(Z- 2.3) 


QO 1.86 1.86 4.2 18.00 16.45 
0.2 fac tal 2.70 4.4 18.90 17-20 
0.4 3044 3.45 4.6 19.80 18.00 
0.6 4.15 4.14 4.8 20.75 18.80 
0.8 4.84 4.80 5) (0 21.75 19.70 
1.0 Seed 0044 5.2 22.75 20.60 
1.2 6.20 6.07 544 23.80 21.50 
1.4 6.88 6.69 Bae 24.90 22 240 
ho Tho TL ToSO 5.8 26.00 23.25 
1.8 slera) 7.95 6.0 27.10 24.20 
2.0 9.04 8.59 6.2 28.20 25.20 
Bow 9.84 9.25 6.4 29.50 26.20 
2.4 10.64 SIRE: 6.6 30.70 27.25 
2.6 11.49 10.60 Ae 32.00 28.30 
208 12.350 11.30 7.0 33.40 29.45 
3.0 13.05 12.00 7.2 34.80 30.55 
30d 13.85 12.795 7.4 26.110 Rall og Ws 
304 14.65 13.50 7.6 37.60 33.00 
3.6 15.45 14.25 7.8 39.10 34.20 
508 16.30 14.95 8.0 40.80 35.50 
4.0 17.15 15.70 
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